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Abstract 

The linear electrodynamic properties of graphene - the frequency-dependent con- 
ductivity, the transmission spectra and cohective excitations - are briefly outlined. 
The non-linear frequency multiplication effects in graphene are studied, taking into 
account the influence of the self-consistent-field effects and of the magnetic field. 
The predicted phenomena can be used for creation of new devices for microwave 
and terahertz optics and electronics. 
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1 Introduction 



Graphene is a new material experimentally discovered about three years ago 
[Tp] . This is a monolayer of carbon atoms packed in a dense two-dimensional 
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(2D) honeycomb lattice, Figure [H left panel. The spectrum of electrons in 
graphene 



Em = (-l)'A 



/3, / = 1,2, 



can be calculated in the tight-binding approximation taking into account the 
symmetry of the lattice (here k = {kx, ky) is the electron wavevector; the 
vectors bj are defined in the caption of Figure [1]). It consists of two bands 
(/ = 1 and 2) which touch each other at the six corners of the hexagon shaped 
Brillouin zone, Figured], right panel. In uniform and undoped graphene at zero 
temperature, the lower band i?ki is fully occupied while the upper band i?k2 
is empty, and the Fermi level goes through these six, so called Dirac points 
Qj, i = 1, . . . , 6. Near the Dirac points the electron dispersion is linear, 

EM = i-iyV\k\, k = k-Qi, |k|a<l, 1/ ~ 10^ m/s, (2) 



and the behavior of electrons can be described by the Dirac equation with the 
Hamiltonian H^, = VaaPa, where o"q, are Pauli matrices {a = x,y). Only two 
of the six Dirac cones are physically inequivalent which is accounted for by 
the valley degeneracy factor = 2. 

The linear (massless) dispersion of graphene electrons near the Fermi level 
leads to a number of interesting linear and non-linear electrodynamic phe- 
nomena, which are briefly discussed in this paper. Some of these effects are 
very promising for microwave and terahertz applications of graphene. 
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2 Linear electromagnetic response 



2.1 Frequency- dependent conductivity 

Using the standard self-consistent linear-response theory and starting from 
the Dirac Hamiltonian Ho one can calculate the frequency dependent con- 
ductivity of graphene criuj) = cTintrai^^) + o'interi^^) [3]fill5]|6ll71l8] ■ It consists of 
two contributions. The intra-hand conductivity 

has the Drude form and is similar to the conductivity of conventional 2D 
electron systems. The mter-band conductivity has the form 

^ _ ie^gsOv 7 sinhx h{uj + i'y)/2T ^ 

ainter[UJ) - J ^^^^ ^ ^^^^ ^ X _ ^^^^ ^ ^^)/2r]2 

where T, fi and 7 are the temperature, the chemical potential and the mo- 
mentum scattering rate respectively, and (7^ = 2 is the spin degeneracy. Figure 
12^ shows the total graphene conductivity as a function of u/fi aX several 
typical values of T//i and 7//X. At high frequencies, u ^ max{T, /i}, ct{uj) 
tends to a universal value cr^^oo e'^/ih dependent only on the fundamen- 
tal physical constants. The behaviour of cr(co') shown in Figure [2^ has been 
experimentally confirmed in Ref. [9]. 

2.2 Transmission, reflection and absorption of light 

The experimentally measurable transmission T = reflection R = |rp and 
absorption A = 1 — T — R spectra are determined by the conductivity o'{u!), 
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since the transmission amplitude tist = l + r=[l + 27rcr(ct;)/c] ^. At high 
frequencies the transmission and absorption coefficients are universal, 

T = l- yl=(l + 7ra/2)-^ ^ 1 - na, u; > max{T, /x}, (5) 

where a = /he. Figure [2t> illustrates the Tiu) dependence [S]. The universal 
transmission ([5]) has been observed in [TO] . 

2.3 Plasmons and other eleetromagnetie exeitations 

The plasma waves in graphene have been studied in Refs. p!T|12|13j . In the 
long- wavelength limit q <^ kp the 2D plasmons in graphene have a square-root 
dispersion {kp is the Fermi wave-vector) 

i^l = gsgve'^fiq/2h^ oc y/rTsq- (6) 

The density (n^) dependence of the 2D plasmon frequency in graphene Up oc 
differs from that of the conventional 2D plasmons {cu oc n]^'^). At larger 
wavevectors {q ~ kp, but still qa <^ 1) the 2D plasmons acquire an additional 
damping due to the inter-band absorption [12j. 

At even larger wavevectors q ~ Qj (g ^ 1/a) the 2D plasmon spectrum 
should be calculated taking into account the full energy spectrum of graphene 
electrons ([T]) and the local field effects [1^ . Such calculations [15] reveal a new 
type of low-frequency plasmons - the inter-valley plasmons - with the linear 
dispersion u = S\q — QJ, where i = 1, . . . ,6 and S >V. 

The 2D plasmons ([6j) are the transverse magnetic (TM) modes. The transverse 
electric (TE) electromagnetic modes do not exist in conventional 2D electron 
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systems. It was shown however [8] that in graphene the TE modes should exist 
at < 2yU, where the imaginary part of the dynamical conductivity (j"{uj) is 
negative, Figure 



3 Non-linear electromagnetic response 



3.1 Frequency multiplication 



If a particle with the linear dispersion ([2]) is placed in the external electric 
field E(t) = Eo cos cut, its momentum p{t) will be proportional to sinut and 
the velocity v(t) - to v = Vp-E'(p) = V^p/|p| ~ sgn(sinci;t). As the function 

, , 4 / 1 1 \ , , 

sgn(sin cut) = — sin cut + - sin Scut + - sin Scut + . . . (7) 
vr V 3 5 / 



contains all odd Fourier harmonics, irradiation of graphene by a wave with the 
frequency uj should lead (in contrast to the conventional electron systems) to 
the emission at higher harmonics at frequencies muj with m = 1, 3, 5, . . . 



A more accurate theory [T7] takes into account the distribution of electrons 
over quantum states in the energy bands and the self-consistent field effects. 
If T <C /i, the higher harmonics generation depends on two parameters S = 
eEoV/ufi and T/u, where F = 2afi/h is the (linear-response) radiative decay 
rate [17] (here a = e^/tic). Our calculations show (Figure [3]) that the system 
efficiently generates higher harmonics ii £ ^ max{l,F/fi}, i.e. at 

So > ^ = — ~ 300— X lO^Vcm^ . 8 

eV c cm 
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The estimate ([8]) does not depend on the frequency of radiation (it is assumed 
however that the dimensions of the sample exceed the wavelength of radiation 
and that hw <2fi). 



3.2 Response to a pulse excitation 

Electromagnetic response of graphene to a strong pulse excitation E(t) = 
Eoro5(t) also differs from that of conventional 2D electron systems (here Eq 
and To are the amplitude and the duration of the pulse). The momentum re- 
laxation in graphene after the strong pulse excitation is linear in time [17] , in 
contrast to the exponential relaxation in conventional systems. The charac- 
teristic response time in the non-linear regime is ~ eE^ToV / . 

3.3 Non-linear response in a magnetic field 

To describe the influence of a magnetic field B = (0, 0, B) on the non-linear 
electromagnetic response of graphene we solve the system of quasi-classical 
non-linear equations of motion for ^ 1 Dirac quasi-particles 



Here Vj = [px + ipy)j/pF is the complex momentum of the j-th particle, 
normalized to the Fermi momentum pp, J = [j^ + ijy)/ensV is the current, 
T = ujt and B = eBV"^ /ufic. The term TJ/uj in ([9]) results from the self- 
consistent-field effects and describes the radiative decay. 

Figure shows the time dependence of the ac electric current in the regime 
of the strong ac electric £ ^ 1 and weak magnetic fields B <^ £. Apart from 



dr 




(9) 
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the current jx{t) with the time dependence close to ([7]) the system generates 
the Hall current jy{t) with the higher frequency harmonics. In the regime of 
weak electric fields S B the frequency transformation effects become even 
more complicated, Figure |Dd. The regular particle dynamics becomes chaotic 
and the system generates both lower and a lot of higher harmonics. This is a 
consequence of the singular motion of graphene electrons in the vicinity of the 
Dirac point where the cyclotron frequency Uc = eBV/pc of individual particles 
diverges at p — 0. The transition to the chaotic dynamics should be observed 
at Eo < VB/c. 



4 Summary 



Due to the linear energy dispersion ([2]), graphene is a strongly non-linear 
material from the viewpoint of its electrodynamic properties. The frequency 
multiplication effect falls down very slowly with the harmonics index m and 
should be seen in moderate electric fields. In weak magnetic fields it should 
be possible to observe the transition from the regular to chaotic particle dy- 
namics dependent on the amplitudes of the B and E fields. The predicted 
non-linear electrodynamic phenomena open up new interesting opportunities 
for studying the fundamental physics of Dirac quasi-particles as well as for 
building innovative devices for microwave and terahertz optoelectronics. 

I thank Igor Goychuk and Timur Tudorovskiy for useful discussions. The work 
was supported by the Swedish Research Council and INTAS. 
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Fig. 1. Left: The honeycomb lattice of graphene. The basis lattice vectors are 
ai = (-v/3/2, — l/2)a and a2 = (0, l)a; the b vectors are hj = {a/2){cos 4>j, sin 

= tt/3 + JItt/S, j = 0,1,2, where a = 2.46A is the lattice constant. Right: 
The band structure of graphene electrons ([1]). The energy is normalized to the 
full width of the energy band A = 2\/^hV/a. The six Dirac points lie at 
Qj = (47r/3a)(cos (pj, sin cpj), (pj = ir/Q + j'vr/S, j = 0, . . . , 5. 
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Fig. 2. (a) The real and imaginary parts of the conductivity (7{uj), in units 
and (b) the transmission coefficient of graphene as a function of fko/fi. 
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eE(|V/a)n r/m 



Fig. 3. The Fourier components of the current j/eUgV as a function of (a) the field 
parameter 6 = eEoV/uiJ, at T/uj = and (b) of F/o; at 6 = 10. 
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Fig. 4. Electric current j{t)/ensV vs time at (a) £ = 5 and (b) £ = 0.1 and at 

;S = l,r/u; = l. The E-Held is polarized in the x-direction. 
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